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Indirect support for the evidence of neutrinoless double beta decay 
C$ • reported recently, is obtained by analysis of other Ge double beta 

experiments, which yield independent information on the background 
in the region of Qpp. Some statistical features as well as background 
simulations with GEANT 4 of the HEIDELBERG-MOSCOW expe- 
riment are discussed which disprove recent criticism. 
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Recently first experimental evidence has been reported for neutrinoless do- 
uble beta decay. Analysis of 55 kg y of data, taken by the HEIDELBERG- 
MOSCOW experiment in the GRAN SASSO over the years 1990 - 2000, has 
led P El H to a half-life 



T 1/2 = (0.8 - 18.3) x 10 25 years (95% C.L.), 



(1) 



with best value of Ti/ 2 = 1.5 x 10 25 y, for the decay of the double beta emitter 
76 Ge 



76 



Gc 



.76 



Se + 2e~. 



(2) 



Assuming the decay amplitude to be dominated by exchange of a massive 
Majorana neutrino (see, e.g. [Hj), this half-life results in a value of the effec- 
tive neutrino mass 



< m >= \H U l 



0.05 - 0.84 eV (95% C.L. 



(3) 



with best value of 0.39 eV. Here a 50% uncertainty in the nuclear matrix 
elements has been taken into account (for details see jHj). 
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Figure 1: The spectrum taken with Ge detectors Nrs. 1,2,3,4,5 over the period 
August 1990 - May 2000 (54.98 kg y), in the energy range 2000 - 2100 keV. Simul- 
taneous fit of the 214 Bi-lines and the two high-energy lines yield a probability for 
a line at 2039 keV of 91% C.L. 
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Figure 2: Result of the peak-searching procedure performed on the 
HEIDELBERG-MOSCOW spectrum (taken with detectors 1,2,3,5) using the Ma- 
ximum Likelihood approach (left) and the Bayesian method (right). On the y axis 
the probability of having a line at the corresponding energy in the spectrum is 
shown. 
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Figure 3: Result of the peak-search procedure performed for the UCBS/LBL 
spectrum |22j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



This is for the first time that the absolute scale of the neutrino mass spectrum 
has been fixed, which cannot be achieved by neutrino oscillation experiments. 
This result restricts possible neutrino mass scenarios to degenerate or (still 
marginally allowed) inverse hierarchy El Ej • In the degenerate case it leads 
to a common neutrino mass eigenvalue of 

m x =0.05 -3.4 eV (95% C.L.). (4) 

This result is nicely consistent with later collected or analyzed experimental 
data, such as results from Large Scale Structure and CMB measurements 
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Figure 4: Analysis of the spectrum measured by D. Caldwell et al. |22| . with 
the Maximum Likelihood Method, in the energy range 2000 - 2060 keV assuming 
lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0 keV. No indication for a signal 
at 2039 keV is observed in this case. 



PIU ITT| IT2"] , or ultra- high cosmic rays J3| • The former yield an upper limit 
of J2i m i=l--0 eV (corresponding in the degenerate case to a common mass 
eigenvalue mo < 0.33 eV). The Z-burst scenario for ultra-high cosmic rays 
requires 0.1 - 1.3 eV [13]. Tritium single beta decay cuts the upper range in 
eq.© down to 2.2 or 2.8 eV [TJj. 

There is further theoretical support for a neutrino mass in the range fixed by 
the HEIDELBERG-MOSCOW experiment. A model basing on an A4 sym- 
metry of the neutrino mass matrix requires the neutrinos to be degenerate 
and the common mass eigenvalue to be >0.2 eV . Starting with the hy- 
pothesis that quark and lepton mixing are identical at or near the GUT scale, 
Mohapatra et. al. ^H] show that the large solar and atmospheric neutrino 
mixing angles can be understood purely as result of renormalization group 
evolution, if neutrino masses are quasi-degenerate (with same CP parity). 
The common Majorana neutrino mass then must be larger than 0.1 eV. The 
fact that WMAP and less stringent tritium decay cuts away the upper part 
of the allowed range from double beta decay (eq. (4)), may indicate [S] that 
indeed the neutrino mass eigenvalues have the same sigen of CP phases - as 
required by (THj. 

In this Letter we report additional support of the double beta result of [U El 
El El 5 by further discussion of the structure of the experimental background, 
by statistical considerations and by analysis of other double beta experiments 
investigating the decay of 76 Ge. 

Important points in the analysis of the measured spectrum are the structure 
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of the background around Qpp (=2039.006(50) keV [Hj), and the energy 
range of analysis around Qpp. 

Background lines in the vicinity of Qpp: Fig. 1 shows the spectrum mea- 
sured in the range 2000 - 2100 keV in its original binning of 0.36 keV. By the 
peak search procedure developped El on basis of the Bayes and Maximum 
Likelihood Methods, exploiting as important input parameters the experi- 
mental knowledge on the shape and width of lines in the spectrum, weak 
lines of 214 Bi have been identified at the energies of 2010.78, 2016.7, 2021.6 
and 2052.94 keV P El El i]- Fig. 2 shows the probability that there is a 
line of correct width and of Gaussian shape at a given energy, assuming all 
the rest of the spectrum as flat background (which is a highly conservative 
assumption). 

The intensities of these lines have been shown to be consistent with other, 
strong Bi lines in the measured spectrum according to the branching ratios 
given in the Table of Isotopes |18j . and to Monte Carlo simulation of the 
experimental setup j^j. Note that the 2016 keV line, as an E0 transition, can 
be seen only by coincident summing of the two successive lines E = 1407.98 
keV and E = 609.316 keV. Its observation proves that the 238 U impurity 
from which it is originating, is located in the Cu cap of the detectors. Recent 
measurements of the spectrum of a 214 Bi source as function of distance 
source-detector confirm this interpretation [TH] . 

Premature estimates of the Bi intensities given in [23] thus are incorrect, 
because this long-known spectroscopic effect of true coincident summing [21] 
has not been taken into account, and also no simulation of the setup has 
been performed (for details see jSlED])- 

These 214 Bi lines occur also in other investigations of double beta decay with 
Ge - and - even more important - also the additional structures in Fig. 2, 
which cannot be attributed at present, are seen in these other investigations. 
There are three other Ge experiments which have looked for double beta 
decay of 76 Ge. First there is the experiment by Caldwell et al. [22], using 
natural Germanium detectors (7.8% abundance of 76 Ge, compared to 86% 
in the HEIDELBERG-MOSCOW experiment). This was the most sensitive 
natural Ge experiment. With their background a factor of 9 higher than 
in the HEIDELBERG-MOSCOW experiment and their measuring time of 
22.6 kg years, they had a statistics of the background by a factor of almost 
four larger than in the HEIDELBERG-MOSCOW experiment. This gives 
useful information on the composition of the background. 
Applying the same method of peak search as used in Fig. [21 yields indica- 
tions for peaks essentially at the same energies as in Fig. [21 (see Fig. EJ). This 
shows that these peaks are not fluctuations. In particular it sees the 2010.78, 
2016.7, 2021.6 and 2052.94 keV 214 Bi lines, but a 1 s o the unattributed lines 
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Figure 5: Result of the peak-search procedure performed for the ITEP/YePI 
spectrum |23j (left: Maximum Likelihood method, right: Bayes method). On the 
y axis the probability of having a line at the corresponding energy in the specrtum 
is shown. 
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Figure 6: Result of the peak-search procedure performed for the IGEX spectrum 
|24l I26j using the ML approach (left) and the Bayesian statistics (right). On the y 
axis the probability of having a line at the corresponding energy in the spectrum 
is shown. 



at higher energies. It finds, however, n o line at 2039 keV. This is consistent 
with the expectation from the rate found in the HEIDELBERG-MOSCOW 
experiment. About 16 observed events in the latter correspond to 0.6 ex- 
pected events in the Caldwell experiment, because of the use of non-enriched 
material and the shorter measuring time. Fit of the Caldwell spectrum al- 
lowing for the 214 Bi lines and a 2039 keV line yields 0.4 events for the latter 
(see Fig. 0J. 

The first experiment using enriched (but not high-purity) Germanium 76 
detectors was that of Kirpichnikov and coworkers |2H]. These authors show 
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only the energy range between 2020 and 2064 keV of their measured spec- 
trum. The peak search procedure finds also here indications of lines around 
2028 keV and 2052 keV (see Fig. EJ), but not any indication of a line at 
2039 keV. This is consistent with the expectation, because for their low 
statistics of 2.95 kgy they would expect here (according to HEIDELBERG- 
MOSCOW) 0.9 counts. 

Another experiment (IGEX) used between 6 and 8.8 kg of enriched 76 Ge, but 
collected since beginning of the experiment in the early nineties till shut- 
down in 1999 only 8.8kgyears of statistics [2HEE]- The authors of (2U EEj 
unfortunately show only the range 2020 to 2060 keV of their measured spec- 
trum in some detail. Fig. |U1 shows the result of our peak scanning of this 
range. Clear indications are seen for the Bi lines at 2021 and 2052 keV, but 
also (as this is the case for the spectrum of [2H], see Fig.5 ) of the unidenti- 
fied structure around 2030 keV. Because of the conservative assumption on 
the background treatment in the scanning procedure (see above) there is no 
chance to see a signal at 2039 keV because of the 'hole' in the background 
of that spectrum (see Fig. 1 in [211 )• With some good will one might see, 
however, an indication of 3 events here, consistent with the expectation of 
the HEIDELBERG-MOSCOW experiment of 2.6 counts. 




Figure 7: Example of a random-generated spectrum with a Poisson distributed 
background with 4.0 events per channel and a Gaussian line centered in channel 
50 (line-width corresponds to a standard-deviation of a = 4.0 channels). The left 
picture shows a spectrum with a line-intensity of 10 events, the right spectrum a 
spectrum with a line-intensity of 100 events. The background is shown dark, the 
events of the line bright. 

The power of the peak search procedure: At this point it may be useful to 
demonstrate the potential of the used peak search procedure. Fig. [7| shows 
a spectrum with Poisson-generated background of 4 events per channel and 
a Gaussian line with width (standard deviation) of 4 channels centered at 
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Figure 8: Results of analysis of random- number generated spectra, using Bayes 
and Maximum Likelihood method (the first one with different prior distributions). 
For every number of events in the simulated line, shown on the x-axis, 1000 random 
generated spectra were evaluated with the five given methods. The analysis on 
the left side was performed with an Poisson distributed background of 0.5 events 
per channel, the background for the spectra on the right side was 4.0 events per 
channel. Every vertical line shows the mean value of the calculated best values 
(thick points) with the la error area. The mean values are in good agreement 
with the expected values (black dashed lines). 

channel 50, with intensity of 10 (left) and 100 (right) events, respectively. Fig. 
IS1 right shows the result of the analysis of spectra of different line intensity 
with the Bayes method (here Bayes 1-4 correspond to different choice of the 
prior distribution: (1) nirf) = 1 (flat), (2) nirj) = I/77, (3) fi(r]) = 1/y/rj, (4) 
Jeffrey's prior) and the Maximum Likelihood Method. For each prior 1000 
spectra have been generated with equal background and equal line intensity 
using random number generators available at CERN |28j . The average values 
of the best values agree (see Fig. |HJ) very well with the known intensities also 
for very low count rates (as in Fig. left). 

A detailed simulation has been done to show the behavior of the position of 
the peak maximum of a line of 15 events, with width (standard deviation) of 4 
channels, on a Poisson-distributed background with 0.5 events/channel. For 
that reason, a peak scan procedure was performed on 1000 randomly created 
spectra. For each spectrum, the deviation of the real position of the line to 
the position of the maximum probability was determined. The distribution 
of these deviations is shown in figure^! which gives a RMS value of 4.00 chan- 
nels, corresponding to ~ 1.44 keV in the spectra from the HEIDELBERG 
MOSCOW experiment. The figure M describes the possible degree of devia- 
tion of the energy of the peak maximum from the transition energy, on the 
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Figure 9: Simulated spectra with a Poisson-distributed background and a Gau- 
sian line with 15 events centered at channel 50 with a width (standard deviation) of 
4.0 channels, created with different random numbers. Shown is the distribution of 
the deviation of the real position of the lines, obtained by the peak scan procedure, 
for 1000 spectra. The rms value of this distribution is 4.00 channels, corresponding 
to 1.44 keV in the spectra of the HEIDELBERG-MOSCOW experiment. 



level of statistics collected in experiments like the HEIDELBERG-MOSCOW 
experiment. This should be considered when comparing Figs. 121315101 
Influence of the choice of the energy range of analysis: The influence of the 
choice of the energy range of the analysis around Qpp has been thoroughly 
discussed in [2J E] ■ Since erroneous ideas about this point are still around, a 
few further comments may be given here. In Fig. El we show the analysis 
of a simulated spectrum consisting of a Gaussian line of width (standard 
deviation) of 4 channels and intensity of 10 counts on a Poisson-distributed 
background of 4 events per channel (see fig. H left) , as function of the width 
of the range of analysis. It is seen that a reliable result is obtained for 
a range of analysis of not smaller than 35 channels (i.e. ±18 channels) - 
one channel corresponding to 0.36 keV in the HEIDELBERG-MOSCOW 
experiment. Fig. ITT1 shows this in a more general way. Every dot represents 
the mean deviation of the calculated line intensity to the real line intensity, 
obtained from 100 random spectra. The exact value for the deviation A is 
calculated by 

A = -jkE|n-&i| AT = 100 n = 0,5,10,20,50, 
i=i 

with fej being the result for the line intensity from the evaluation and n the 
real line intensity, chosen as 10 events in Fig. ^2 The spectra used here have 
again a background of 4.0 events per channel. Also in this general case it is 
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Figure 10: Result of an analysis as function of the evaluation width. The used 
spectrum consists of a Poisson distributed background with 4 events per channel, 
and a line of 10 events (see fig. left part). The dark area corresponds to a 68.3% 
confidence area with the dark line being the best value. Below an evaluation width 
of 35 channels the result becomes unreliable, above 35 channels the result is stable. 

seen that, if the evaluation width is more than 35 channels, the deviation is 
constant, so that an evaluation width of at least 35 channels is required for 
reliable results. This is an important result, since it is of course important to 
keep the range of analysis as s m a 1 1 as possible, to avoid to include lines in 
the vicinity of the weak signal into the background. This unavoidably occurs 
when e.g. proceeding as suggested in j2Hj. The arguments given in those 
papers are therefore incorrect. Also Kirpichnikov, who states [29J that his 
analysis finds a 2039 keV signal in the HEIDELBERG-MOSCOW spectrum 
on a 4 sigma confidence level (as we also see it, when using the Feldman- 
Cousins method |27j). makes this mistake when analyzing the pulse-shape 
spectrum. 

The background around Qpp from GEANT4 simulation: Finally the back- 
ground around Qpp will be discussed from the side of simulation. A very 
careful new simulation of the different components of radioactive background 
in the HEIDELBERG-MOSCOW experiment has been performed recently 
by a new Monte Carlo program based on GEANT4 |3lH l31j. This simulation 
uses a new event generator for simulation of radioactive decays basing on 
ENSDF-data and describes the decay of arbitrary radioactive isotopes in- 
cluding alpha, beta and gamma emission as well as conversion electrons and 
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Figure 11: Simulation of spectra with a line intensity of 10 events for a Poisson- 
distributed background with 4.0 events/channel. The horizontal axis shows the 
energy range of analysis. The vertical axis shows the mean deviation of the inten- 
sities obtained by the peak scanning procedure from the known real intensities. 
Each dot represents the mean deviation obtained from 100 random spectra. The 
deviation is constant, when more than 35 channels are used for the evaluation. 




X-ray emission. Also included in the simulation is the influence of neutrons 
in the energy range from thermal to high energies up to 100 MeV on the 
measured spectrum. Elastic and inelastic reactions, and capture have been 
taken into account, and the corresponding production of radioactive isotopes 
in the setup. The neutron fluxes and energy distributions were taken from 
published measurements performed in the Gran Sasso. Also simulated was 
the cosmic muon flux measured in the Gran Sasso, on the measured spec- 
trum. To give a feeling for the quality of the simulation, Fig. El shows 
the simulated and the measured spectra for a 228 Th source spectrum for 4 
of our five detectors. It should be noted, that this simulation is not a fit. 
It is a calculation of the expected spectra on the basis of the knowledge of 
the geometry of the setup, and its properties of absorption and scattering of 
gamma radiation, and of the known position and strength of the source. The 
agreement is excellent. 

The simulation of the background of the experiment reproduces a 1 1 ob- 
served lines in the energy range between threshold (around 100 keV) and 
2020 keV [23 EH1 ED] • Fig. H31 shows the simulated background in the range 
2000-2100 keV with all k n o w n background components. The black solid 
line corresponds to the measured data in the period 20.11.1995 - 16.4.2002 
(55.57 kg y). It should be noted here, that the simulated spectrum is not 
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Figure 12: Comparison of the measured data (black line, November 1995 to April 
2002) and simulated spectrum (red line) for the detectors Nrs. 1,2,3 and 5 for a 
232 Th source spectrum. The agreement of simulation and measurement is excellent. 



folded with the energy resolution of the detectors (which explains the sharp 
'lines' in the simulated spectrum). 

The background around Qpp is according to the simulations flat, the 
only expected lines come from 214 Bi (from the 238 U natural decay chain) at 
2010.89, 2016.7, 2021.6, 2052.94, 2085.1 and 2089.7 keV. Lines from cosmo- 
genically produced 56 Co (at 2034.76 keV and 2041.16 keV), half-life 77.3 days, 
are not expected since the first 200 days of measurement of each detector are 
not used in the data analysis. Also the potential contribution from decays 
of 77 Ge, 66 Ga, or 228 Ac, should not lead to signals visible in our measured 
spectrum near the signal at Qpp. For details we refer to [3*T] . 
The structures around 2028 keV, 2066 keV and 2075 keV seen - as also the 
214 Bi lines - in practically all Ge experiments (see above), cannot be identified 
at present. 

Conclusion: Concluding, additional support has been given for the evi- 
dence of a signal for neutrinoless double beta decay, by showing consistency 
of the result - for the signal, and for the background - with other double 
beta decay experiments using non-enriched or enriched Germanium detec- 
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Figure 13: Simulated background of the HEIDELBERG-MOSCOW experiment 
in the energt range from 2000 to 2100 keV with all known background components. 
The black trough-drawn line corresponds to the measured data from 20.11.1995 to 
16.4.2002 (49.59 kg y). 



tors. In particular it has been shown that the lines seen in the vicinity of 
the signal (including those which at present cannot be attributed) are seen 
also in the other experiments. This is important for the correct treatment of 
the background. Furthermore, the sensitivity of the peak identification pro- 
cedures has been demonstrated by extensive statistical simulations. It has 
been further shown by new extensive simulations of the expected background 
by GEANT4, that the background around Qpp should be flat, and that no 
known gamma line is expected at the energy of Qpp. The 2039 keV signal is 
seen o n 1 y in the HEIDELBERG-MOSCOW experiment, which has a factor 
of at least 10, but in general much more, statistics than all other double 
beta experiments. 
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